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1. Introduction 

A many-body system rapidly brought out of an initial equilibrium state by quenching it 
to either a critical point or into a coexistence region of the phase-diagramme where 
there are at least two equivalent equilibrium states undergoes ageing, that is the 
physical state depends on the time since the quench was performed and hence time- 
translation invariance is broken. Furthermore, either the competition between the 
distinct equilibrium states or else non-equilibrium critical dynamics leads to a slow 
dynamics, where observables depend non-exponentially on time. Both aspects of ageing 
can be conveniently studied through the two-time response and correlation functions 
defined as 

0) . /«(.)-^».- (1.1) 

Cit,s) = m,r)<f^{s,r)) = s-'fc(^^ , /c(y) y"'^/^ (1.2) 

where r) is the order-parameter at time t and location r and h is the conjugate 
magnetic field. These dynamical scaling forms ( ll.l|ll.2p are expected to hold in the 
scaling hmit where both t,s ^ tmicro and also t — s ^ ^micro, where tmicro is some 
microscopic reference time. In writing those, it is implicitly assumed that the underlying 
dynamics is such that there is a single relevant length-scale L = L(t) ~ t^^^, where z 
is the dynamical critical exponent. Universality classes of ageing are defined through 
the values of the non-equilibrium exponents a, b, Ac, ^r- for a non-vanishing initial 
value of the order-parameter, it is known that Xc,r may be related to the slip exponent 
introduced in [1] but there are no scaling relations with equilibrium exponents. For 
reviews, see e.g. [21 [3l HI [5], E] . 

A lot of attention has been devoted to the case where the underlying dynamics 
is such that there are no macroscopic conservation laws and in particular the order- 
parameter is non-conserved. In this situation, which is appropriate for the description 
of the ageing of magnets, results for the exponents a, b, \c, Xr have been derived. 
Furthermore, it appears that simple dynamical scaling can be extended to a local scale- 
invariance, where the response functions transform covariantly under more general time- 
dependent rescalings of time and space [8] . In particular, this leads to exphcit predictions 
for the scaling functions (and also for fc{y) if 2; = 2) which have been tested and 

confirmed in a large variety of systems, see [7j for a recent review. 

In this paper, we are interested in carrying out a first case-study on the ageing and 
its local scale-invariance when the order-parameter is conserved. Physically, this may 
describe the phase-separation of binary alloys, where the order-parameter is given by 
the concentration difference of the two kinds of atoms. Another example are kinetic 
growth-processes at surfaces when mass conservation holds. One of the most simple 
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models of this kind is the Mulhns- Herring model (U [TU], which leads to a dynamical 
exponent z = 4 and has recently been discussed in detail in fH] . 

Presently, there exist very few studies on ageing with a conserved order-parameter. 
In a numerical study of the 2D Ising model quenched to T < Tc, where z = 3, 
Godreche, Krzakala and Ricci-Tersenghi [12] calculated the scaling function fc{y) which 
remarkably shows a cross-over between a first power- law decay fc{y) ~ y~^c/^ for 
intermediate values of y to the final asymptotic behaviour (11.21) for y ^ 1. They find a 
value of X'q = 2.5, whereas the value for Ac is much smaller with Ac = 1- On the other 
hand. Sire [13] calculated fc{y) exactly in the critical spherical model (T = Tc), where 
z = 4, and found a single power-law regime|il For a review, see [4j. 

Here, we shall revisit the spherical model with a conserved order-parameter, 
quenched to T = Tc. Besides reviewing the calculation of the two-time correlation 
function we shall also provide the exact solution for the two-time response function. This 
will be described in section 2. Our main focus will be to show how these exact results can 
be understood from an extension of dynamical scaling to local scale-invariance (LSI) with 
z = 4. In section 3, we recall the main results of LSI for that case and in particular write 
down a linear partial differential equation of which this LSI is a dynamical symmetry. 
In section 4, we consider the Langevin equation which describes the kinetics with a 
conserved order-parameter (model B in the Halperin-Hohenberg classification) and show 
for the case at hand how the Langevin-equation can be split into a 'deterministic' and a 
'noise' part, so that all averages to be calculated can be exactly reduced to quantities to 
be found in a noise-less theory. Local scale-invariance directly determines the latter, as 
we shall see in section 5, where we can finally compare with the exact spherical model 
results of section 2. This is, together with the earlier study of the MuUins-Herring 
model in [11], the first example where local scale-invariance can be confirmed for an 
exactly solvable model with dynamical exponent z ^ 2. Our conclusions are formulated 
in section 6. Some technical points are reported in the appendices. In appendix A 
we derive the two-point function from the symmetry operators and in appendix B we 
calculate solution of an integral which appears frequently in our computations. In an 
outlook beyond the present specific context, we discuss in appendix C the autoresponse 
function for values of z different from 4. In appendix D we discuss the most simple 
possible model, namely a free particle submitted to conserved noise. This serves as 
a further illustration of some of the space-time properties of the conserved spherical 
model. 

I It is well-known that the spherical model with conserved order-parameter quenched to below Tc 
shows a multiscaling which is not captured by the simple scaling form ()1.2p [Mj . It is understood that 
this is a peculiarity of the n — > oo limit of the conserved 0(n) vector-model [15) . 
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2. The spherical model with a conserved order-parameter 

The spherical model was conceived in 1953 by Berlin and Kac as a mathematical model 
for strongly interacting spins which is yet easily solvable and it has indeed served a 
useful role in providing exact results in a large variety of interesting physical situations. 
It may be defined in terms of a real spin variable S{t, x) attached to each site of a 
hypercubic lattice A <Z Z'^ and depending on time t, subject to the (mean) spherical 
constraint 

/j2S{t,xA=Ar (2.1) 

XccGA / 

where M is the number of sites of the lattice. The Hamiltonian is 7-^ = — J2{xy) ^x^y 
where the sum is over pairs of nearest neighbours. The kinetics is assumed to be given 
by a Langevin equation of model B type 

dtS{t, x) = -Vl [Visit, x) + 3(t)S(t, x) + h{t, x)] + r]{t, x) (2.2) 

where ^{t) is the Lagrange multiplier fixed by the mean spherical constraintl§| and the 
coupling to the heat bath with the critical temperature Tc is described by a Gaussian 
noise rj of vanishing average and variance 

{r]{t, x)r]{t', x')) = -2T,Vl6{t - t')5{x - x'). (2.3) 

and h(t, x) is a small external magnetic fieldjjj] Here the derivative on the right-hand side 
of (12. 3p expresses the fact that the noise is chosen in such a way that it does not break 
the conservation law. On the other hand, there are physical situations such as surface 
growth with mass conservation, where this is not the case. The non-conserved MuUins- 
Herring model [9l [10] is given by equation (12. 2p with 3 = and the noise correlator 
{r]{t,x)ri{t',x')) = 2Tc6{t — t')5{x — x'), see [T^ for details. Here we shall study the 
conserved Mullins-Herring model, that is equation (12. 2p with ^{t) = and the noise 
correlator (12. 3p . A similar remark applies on the way we included the perturbation 
h{t, x). Here we study the case, where the perturbation does not break the conservation 
law. 

The spherical constraint is taken into account through the Lagrange multiplier 
3(t) which has to be computed self-consistently. These equations have already been 
considered several times in the literature [HI [IHl |13] . 

§ Considering the spherical constraint in the mean simphfies the calculations and should not affect the 
scaling behaviour. See I16j for a careful study of this point in the non-conserved case. 
II Since we are only interested in averages of local quantities and the initial magnetisation is assumed 
to vanish, this description is sufficient and we need not consider the analogue of the more elaborate 
equations studied by Annibale and Sollich [17j for the non-conserved case. 
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We repeat here the main steps of the calculation of the correlation function [T9l [T3] . 
Equation (12.21) with h{t, a;) = is readily solved in Fourier-space, yielding 



S{t,k) 



-1 

S{0,k)+ J dr exp(u{T,k)^ri{T,k) exp - fc) j (2.4) 
where S{t, k) denotes the spatial Fourier transformation of S{t, r) and 

Lo{t,k) = kH-k^ [ dr^ir), with k := \k\ (2.5) 



Next, uj{t, k) is determined from the spherical constraint. This is easiest in 
the continuum limit where spatial translation-invariance leads to {S^it,x)) = 
(27r)-'^/ dfe(^(t,-fc)^(t,fc)) = 1 where A is the inverse of a lattice cutoff. This has 
already been analysed in the literature [T8|, [191 IE], with the result 

u}{t, k) = kH - Qdkh^'^ (2.6) 

where the constant is determined by the condition 

(2.7) 



^ dxx'^^i |2x2 ^ dy exp [-2x^(1 - y) + 2Qdx\l - ^)\ - ^} = 



In particular, fld = for d > 4 and asymptotic forms for d 2 and d ^ A are 
listed in [18]. Then from (12. 4p the two-time correlator C{t,s,k) := (^S{t, —k)S{s, k) 
is straightforwardly derived 



(2.8) 



C{t, s, k) = so(fc)^ exp ^ — uit, k) — uo{s, k) 
+ 2Tc J du k"^ exp ^ — a;(t, k) — ijj{s, k) + 2uj{u, fc) j 

where so(fc)^ := <y5'(0, — fc)S'(0, k)^. Transforming back to direct space, one obtains 

C(t, s- r) = Ai(t, s- r) + A2(t, s- r) (2.9) 

where v4i(t, s; r) and A2{t, s\ r) are the inverse Fourier transforms of the first and the 
second line of the right-hand side of equation (12. 8p . respectively. They are explicitly 
given by the following expressions, where we use the shorthand y := t/s 

AAt,.,r) = (i + ,)-/*|^^^,„W^exp(-^^) (2.10) 

4 , . ,2 + ! 



X exp —k + Qdk 
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and 



A,{t, s; r) = 2T,.-('^-2)/^ [\e {y + 1 - 26^'^'^' f T^k' {2.11] 

Jo J^d i^^r 



ikr \ f ,^ 2/7/1/2 + 1-201/2 



The term A2{t,s;r) is the dominating one in the scahng regime where t,s ^ oo and 
y = t/s is kept fixed. There we recover the scahng form (11.21) with b = (rf — 2)/4. The 
exphcit expression for the scahng function fc{y) is quite cumbersome but, if the last 
term in the second exponential in eq. (12.111) can be treated as a constant, one has, up 
to normalisation 



This is exact for the spherical model for d > 4, since then = 0. Equation (12.121) 
also applies to the MuUins-Herring equation for any d > 1 (since here 3 = 0, in this 
case 0d = exactly). For y sufficiently large, the above condition is also approximately 
satisfied. In any case, in the limit y ^ 00 we can read off the autocorrelation exponent 
\c = d + 2. 

The scaling behaviour of the spatio-temporal two-time correlator is illustrated in 
figure [H We see that for relatively small values ofy = t/s, the correlation function does 
not decay monotonously towards zero, but rather displays oscillations whose amplitude 
decays with |r|s~i/'i. This behaviour is quite close to the well-known one for the equal- 
time correlation function. As y increases, these oscillations become less pronounced. 
There is no apparent qualitative difference between the cases d < 4 and d > A. 

2.2. The response function 

The response function can be computed in a way similar to the non-conserved case |20j . 
Equation (12.21) is solved in Fourier space yielding an equation similar to (12. 4p 

S{t, k) = S{0, k) + j^dr exp (u{t, k) 

X (^r]{T,k) + k'^h{T,k)^ exp{-uj{t, k)) (2.13) 

From this the Fourier transform of the response function is computed as R{t, s, k) = 
6{S{t, k))/6h{s, —k). Transforming the result back to direct space, we obtain for t > s 

R{t, s; r) = {t- s)-('^+2)/4 f dfc ^2 ( ■ r 



{2ttY ^ V {t-sy/^ 



X exp ( -k^ + 5dk \^'^\,^,, ) (2.14) 
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Figure 1. Scaling of the spatio-temporal correlation function C{t,s;r) ~ A2{t,s;r) 
with s = 20 and for several values of y = t/s for (a) d = 3 (left panel) and (b) = 5 
(right panel). 



It is easy to see that the scahng form fll.2p holds with a = {d — 2)/A = b. For large 
values of y the scaling function fn^y) is given by 



(2.15) 



Again, for d > 4 this expression is exact up to a prefactor. The value of the autoresponse 
exponent therefore is, for d > 2 



X 



R 



d + 2 = X. 



c 



(2.16) 



in agreement with field-theoretical expectations for the 0(n) model 

It is instructive to consider the full space-time response function explicitly for the 
case that %^ = 0, that is for > 4 in the spherical model and for any d in the Mullins- 
Herring equation. The integral in (12.141) is done in appendix B and we find 

^4 



R{t,s;r) 



23d/2^d/2r(^/4) 



{t 



8T 



0-^2 



-(d+2)/4 



3d 1 

+ 



F ( ^ ^ 



2' 4' 256 (t-s) 

4 



2' 4 2' 256 (t-s) 



(2.17) 
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Figure 2. Scaling of tlie spatio-temporal response function R{t, s; r) for s = 20 and 
several values oi y = t/s ioi- d = 3 (left panel) and d = 5 (right panel). 



This expression will be useful later for direct comparison with the results obtained from 
the theory of local scale-invariance. For fl^ 7^ it is more convenient to use the integral 
representation (12.141) . In figure [2] we display the spatio-temporal response function. 
In contrast to the non-conserved case for a fixed value of y = t/s there are decaying 
oscillations with the scaling variable |r|s~^/^. These oscillations disappear rapidly when 
y is increased. This qualitative behaviour also arises, as we show in appendix D, for the 
simple random walk with a conserved noise. 

Having found both correlation and response functions, we can also calculate the 
fluctuation-dissipation ratio X(t,s) = TR(t, s){dC(t, s)/ds)~^ which measures the 
distance from equilibrium [21]. In particular, we find for the fluctuation-dissipation 
limit 



= lim ( lim X{t, s 



oo \t—>oo 



(2.18) 



in agreement with the known result [4j of the conserved n ^ oo limit of the 0(n) model. 
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3. Symmetries of the deterministic equation 



Before we can examine dynamical symmetries in Langevin equations, we need some 
background material on the dynamical symmetries of the deterministic part of such 
equations. Therefore, we recall in this section first some facts about Schrodinger- 
invariance which describes an extension of dynamical scaling for a dynamical exponent 
z = 2. As shown earlier [8], an analogous extension can be built for z ^ 2 and we shall 
write down the relevant generators for this local scale-invariance for z = 4. 

The physical starting point is the covariance of n-point functions such as the 
correlator . . . , ri, . . . , r„, ) = ri) . . . , 0n(^n, rn)) under the rescaling 

of spatial and temporal coordinates 

F^^\b%, bH^; br,, . . . , 6r„) = . . . ri, . . . , (3.1) 

where the Xi are the scaling dimensions of the fields (pi and z is the dynamical exponent. 
By analogy with conformal invariance at equilibrium, one asks whether an extension of 
dynamical scaling to local scale-transformations with a space-time dependent rescaling 
factor b{t, r) might be possible. The most simple case is the one of Schrddinger- 
invariance which describes the dynamical symmetries of the linear diffusion equation 
(or free Schrodinger equation) 

^sch0(r',t) = 0, with S,,^:=\dt-U/l (3.2) 

where the parameter A is related to the inverse diffusion constant. The symmetry group 
of this equation, that is the largest group of transformations carrying solutions of (13.21) 
to other solutions, was already found 1882 by Lie and is the well-known Schrodinger 
group Sch{d), see [22]. It transforms coordinates as {t,r) ^-^ {t',r') = g{t,r) with 

, TZr + vt + a at + p ^~ ~ 
r = ^ , t = ^ ; a6- (3^ = 1 (3.3) 

-ft + 6 -ft + 6 

where 7^ is a rotation matrix and v,a,a,l3,-f and 6 are parameters. Evidently, the 
dynamical exponent is z = 2. The solutions of (13. 2p transform as 

0(t,r) ^ (T,0)(t,r) = fg[g-\t,r)]<P{g-\t,r)) (3.4) 

where the companion function fg is known explicitly [22j. A field transforming as ( 13. 4p 
is called quasiprimary. The generators of the Schrodinger group form a Lie algebra 
scf)^ = \jie{Sch{d)) and it can be shown that a n-point correlator built from quasiprimary 
fields must satisfy the conditions 



n 

E 

i=l 



Xi {<Pi{ti,ri) . . .(j)n{tn,rn)) = Q (3.5) 
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where Xi stands for any generator of 5cf)^ which acts on the field 0j. For example, 
eq. fl3.5p fixes the two-point function of quasiprimary fields completely |23j . 

Following |8], we now write the generalisation of this to the case when z = 4 which 
has already been started in pi]. Consider the 'Schrodinger operator' 

S--=-^9t + ^{VlY. (3.6) 
lb 

which we shall encounter in the context of the conserved spherical model and where 
obviously 2 = 4. Using the shorthands r ■ dr := X]fc=i '^fc f^rfc,V^ := ^^=1 c^r^ 

:= 'Y2k=i '"fc' infinitesimal generators of local scale-transformations read (using a 
notation analogous to P]) 

X_i := -dt 

1 X 

Xq ■= -tdt - - r ■ dr - - 

X, := - t'dt - |t - Ar^ (V^)-i -\tr.dr 

+ 47(r-9.) [vlY' + 2^{d-A){Vl)-' (3.7) 
.— ^.Q^^ _ Tjdn ; where I < i < j < d 

y(*) - -d 
^-1/4 ~ 

Yi% = -t9,-4Ar,(V^)-^ + 87 9,(V^)-' 

where x is the scaling dimension of the fields on which these generators act and 7, A are 
further field-dependent parameters. Here, the generators X-i-i correspond to projective 

(i) ' 

changes in the time t, the generators Y^_-^^^^ are space-translations, generalised Galilei- 
transformations and so on and are spatial rotations. In writing these generators, 
the following properties of the derivative 9" are assumed 

aX^ = a°+/5 , [d^,r] = ad^-^ (3.8) 

and which can be justified in terms of fractional derivatives as shown in detail in 
appendix A of [8]. The operator (V^)~^ can then be defined formally as follows. For 
example, in c? = 2 dimensions, we have 

00 

(V^)-^ := + dlr' ■■= (3-9) 

n=0 

The remaining negative powers of Vr- are then defined by concatenation, e.g. (V^)~^ = 
(V^)~^ ■ (V^)~^. We easily verify the following commutation relations for n G Z 

[(V^)^r,] =n9,(V^)"-^ (3.10) 
[(V^)^ r^] =2n{r- dr) (V^)""^ + n{n-2 + rf)(V^)"-^ (3.11) 



The generators eq. (13. 7p describe dynamical symmetries of the 'Schrodinger 
operator' (13.61) . This can be seen from the commutators of 5* and the generators X 
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from eq. (13 .71) . Straightforward, but a little tedious calculations give, analogously to [S] 

[s, = , [5, r «/j = , [s, yW] = o 

[S,Xq] = - S , [S,R^''^'^] = (3.12) 

This means that for a solution of the 'Schrodinger equation' S(j) = the transformed 
function is again solution of the 'Schrodinger equation'. For the commutator with 
Xi we find 

[5, X,] = -2tS + ^ - ^ - 1 + y ) (3-13) 
hence a dynamical symmetry is found if the field has the scaling dimension 

. = ^ + l-H2 (3.14) 



which for d = 1 reproduces the result of [8]|5[| Generalising from conformal or 
Schrodinger-invariance, quasiprimary fields transform covariantly and their n-point 
functions will again satisfy eq. (13.51) . A quasiprimary field is now characterised by 
its scaling dimension Xi and the further parameters ji, A,. For example, any two-point 
function F^"^^ = (0i02) built from two quasiprimary fields 0i,2 is completely fixed by 
solving the conditions (13. 5p for the generators in (13. 7p . We quote the result and refer to 
appendix A for the details of the calculation. 

F(2)(t _s,x-y):= (0i(t, x)<f)2{s, y)) 

= 5.,.,5ai,-a.5,,-,. {t - s)-^/2 5^c,0(^) {\x-y\{t- s)-^") (3.15) 

Here 0*^'^^ {u) are scaling functions, the Cg are free parameters and the set £' is defined as 
follows 

' {2,4} if ci>4 

^' := <j {2,4,4-4 if 2<rf<4 (3.16) 

{2,4,2-^,4-4 if ^^<2 

We shall see later that boundary conditions may impose further conditions on the Cg. 
The functions (j)^^\u) are given by the series, convergent for all |u| < oo 

oo 



^=0 



with the coefficients 6^ 



^ o4 r(f + i)r(f + i) _ r(£+^-|-i) 



^ For a free field-theory, where x = d/2, this implies 7/A = ^ 
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4.1. General remarks 

Consider the following stochastic Langevin equation 



5*0 = -Y^V^( - + + T] (4.1) 

where the noise correlator respects the global conservation law 

ivir, tMr', t')) = -^V^5(r - r')(5(t - t') (4.2) 

Here and in what follows, we shall often suppress the arguments of the fields for the 
sake of simphcity, if no ambiguity arises. We shall adopt the standard field-theoretical 
setup for the description of Langevin equations, see e.g. pH [25], [26] for introductions. 
The Janssen- de Domincis action can be written in terms of the order-parameter field 
(p and its conjugate response field and reads 



+ dudR^iu, R) (V^0)(m, R) 



(4.3) 



to which an extra term describing the initial noise must be added, by analogy with the 
non-conserved case [27|, [28]H 

J^nu[<f)A] = \j dRdR'^{0,R)(^<f){0,R)<P{0,R)'j^{0,R') (4.4) 

Averages of an observable O are defined as usual by functional integrals with weight 
exp{-J[(f),(j)]), viz. 

{O) := J Vm^]Oexp{~J[<P,^]) (4.5) 

We decompose the action, in the same way as done in [28] for the non-conserved case, 
into a deterministic and a noise part, that is 

^[0,0]= Jo[0,0]+Jb[0] (4.6) 

with 



Jo[0,0] = JdudR ^(^du-^V%{v%-viu))^<p 



(4.7) 



+ It has been shown by Janssen that at the initial time t = 0, the order-parameter field (/)(0, r) and the 
response field (j){0,r) are proportional [53]. 
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and Jb = Jth + Jinit where 

M^] = ^ / d^dR ^{u, R) (V^0) {u, R) (4.8) 

The point of this spht-up is, as we shall show, that the action J7o[0, 0] has nontrivial 
symmetry properties, in contrast to the full action 0], where these symmetries are 
destroyed by the noise. We call the theory with respect to J7o[0, 4>] noise-free and denote 
averages taken with respect to J7o[0,0] only by (. . .)o. Averages of the full theory can 
then be computed by formally expanding around the noise-free theory 

o) = (Oexp(-Jb[0]))^ (4.9) 

The noise-free theory has a Gaussian structure, if we consider the two-component field 
\I' = . Then the factor exp(— J7o[0, 0]) can be written as exp(— f dudRdu'dR' \E'*A^) 
with an antidiagonal matrix A. From this, two important facts can be deduced: 

(i) Wick's Theorem holds [30]. That is, we can write the 2n-point function as 

'0(tl,ri)...02n(t2n,r'2n))^ = ^ (4.10) 

all possible pairings 
P of 1,2, ...,2n 

(pPi{tp„rp^)^P^{tp^,rp^)'^ ^. . . {(t)P2,..Aip2n-i^r p^^_;)(pp^^{tp^^,r p^S) ^ 

[ii) We have the statement 

= (4.11) 



unless n = m. This is due to the antidiagonal structure of A and can be seen 
by performing the Gaussian integral (see for instance [26j, chapter 4) explicitly 
and taking care of the fact, that the inverse matrix A~^ is antidiagonal again. 
Formally, eq. (14.111) is analogous to the Bargman superselection rule which is used 
in Schrodinger- invariant theories with z = 2 [28] . 

With these tools at hand we cam demonstrate, quite analogous to the non-conserved case 
[28] an exact reduction of any average to an average computed only with the noiseless 
theory. For example, for the two-time response function, we have 

R{t,s) = (nms)) = (0(t)0(s)e-'^^[^l)^ = (0(i)0(s))^ (4.12) 

In going to the last line, we have expanded the exponential and use that because of 
Wick's theorem any 2n-point function can rewritten as a sum over a product of n two- 
time functions which in turn are determined by the Bargman rules. From the special 
structure of Jb[<P] it follows that only a single term of the entire series remains. As a 
consequence, the two-time response function does not depend explicitly on the noise. A 
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similar result holds for the correlation function, see section 5. We shall now first consider 
the noise-free theory and find the two-point function from the dynamical symmetry. 
Afterwards, we shall show that the two-point functions of the full noisy theory can be 
reconstructed from this case. 



4-2. The response function of the noise-free theory 



First, we consider the linear response function of the order-parameter with respect to 
an external magnetic field h 



6h{s,y) 



(4.13) 

h=0 

As already mentioned in the previous section, we choose here, and in contrast to |llj . 
a perturbation respecting the conservation law, which means that we have to add the 
term V|j/i to the Langevin equation (14.11) or the term J dRdu hV^cf) to the action (14. 3p . 
Then it is easy to see that the response function (I4.13p is given by 



R(t, s;x-y) = ^</>(t, a^)VXs, y)^ (4.14) 

We make the important assumption that the field (j){t,r), characterised by the 
parameters A and 7 and the scaling dimension x, is quasiprimary. As suggested in 
[HI [28] we consider also the response field (j)(t, r) to be quasiprimary, with parameters 
A = —A and 7 = —7 but the same scaling dimension x = x. 

We can concentrate on the model with f = for the following reason: Suppose 
4>{r,t) is a solution of equation (14. ip with rj = 0. Then we define 

^{t, r) := exp ^ d™(r)V^^ 0(t, r). (4.15) 

\E'(t, r) fulfils the following dynamical equation 

dt^{t,r) = ^Vt.^{t,r), (4.16) 

which is the same equation with v = 0. If suffices thus to consider the problem with 
V = and then to apply the inverse of the gauge transformation (I4.15P for treating the 
case V ^ 0. In this way, the breaking of time-translation invariance is implemented. 



4.3. The case v = 



This case is relevant for d > 4 in the spherical model and for any d in the Mullins- Herring 
model. We compute (I4.14p for the noise-free theory, taking translation-invariance into 
account. Then the response function is given by 



Ro{t,s; r) = VlF^^\t,s,r) 



(4.17) 
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where the two-point function F^'^\t, s,r) has been computed in the last section. We 
obtain, with the scahng variable u = \r\(t — s)~^/^ 

= {t- J2 Cs4>^'\u) (4.18) 

se£" 

where are constants, S" the set of admissible values for s given below, and the solutions 
ai-g given by 



oo 



(4-19) 



with the coefficients 



A priori, s could take the values s G {2, 4, 2 — d, 4 — d} \ {—2, —4, . . .} as derived in 
appendix B. However, certain values of s have to be excluded, since the solution has to 
be regular for m — and has to vanish for u oo. By inspection of the coefficients b\ 
one finds, in a similar way as done at the end of appendix B: 

{2,4} if rf>4 

£" ■= { {2,4,2-rf} if 2<d<4 (4.21) 

{2,4,2 -rf,4- 4 if d<2 

In the specific example of the spherical model in turns out that ^ = |- Then the solution 
with s = A — d & £" disappears and is no longer admissible for d < 2. 

For completeness, we rewrite the solutions 0^*-'(-u) as hypergeometric functions for 
the admissible s, suppressing some constant prefactors which can be absorbed into the 
constants c^. 

r V ; ^ ^ V ^ 4 A' 2' 4' 2 4' 16 y ^ ^ 

X / \u'^\^^'^ ^/S d 731 d ^ d Am^A 

0(4) = i? _ + _ _ I; _ + _ 1 + 4.23 

^ ^ ' V 16 y ^ V2 4 A 2' 2 4' 4' 16 7 ^ ' 

d^'-'\u)= (-^y~''\Fj2-l;^--^,2-^^,~^) (4.24) 
^ ^ ' V 16 y ^ V A 2 4' 4' 2' 16 y ^ ' 

J<-.(„)= f_^V""'\fJi_2;l_^.!_f?,l;_^) (4.25) 

^ ^ ' V 16 / ^ V2 A' 4' 2 4' 2' 16 y ^ ^ 

The constants are not completely arbitrary for the case A < 0, but have to be arranged 
so that 4>{u) — i> for M — > oo. From [311 132], one knows the asymptotic behaviour of 
the hypergeometric functions. In general, there is an infinite series of terms growing 
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exponentially with u, together with terms falling off algebraically. The leading term of 
the exponentially growing series can be eliminated by imposing the following condition 
on the coefficients 



^ r(i)r(f)r(i + f) r(|)r(i + f)r(i + f) 



^ r(|-f)r(2-f)r(|) , r(i - |)r(| - |)r(i) _ ^ 

Here, some of the constants might have to be set to zero, if the corresponding value 
of s is not admissible. Indeed, the condition (14.260 is sufficient to cancel the entire 
exponentially growing series and the remaining part decreases algebraically, but we 
shall not prove this here. The most important case for us is C2~d = c^-d = 0. In this 

case f l4.26p implies the relation C2 = "fTdT^TT^ ^"^^ ^^^y show that (see 
appendix B for the details) 

/4j((,s;r) = V5Fl2'((,s,r) 

= r„(( /^(fc^)=-^exp(-i^^A_L_)exp(-fcT (4.27) 

This prediction of LSI with 2; = 4 is perfectly consistent with the exact results (12.141) 
and (12.171) of the conserved spherical model for the case (i > 4 if we set ^ = | and 
X = |. The case c? < 4 in the spherical model (where we have w 7^ 0) will be treated in 
the next subsection. 



4.4- The case w 7^ 

In this case, the response function is given by 

R{t,s-r) = eMit,sVl)^lF^^\t,s,r) (4.28) 

where we have defined S,t,s '■= /s*'^'^^('^) using the fact that (p is characterised by 
the parameters —A and —7. Therefore, remembering u = \r\(t — s)~^^^ 

Roit, s, r) = ro{t - s)-(^+^)/^ ^^^^^^^^y^(V^)" c.0W(«) (4.29) 

n=0 ■ se£" 

where = u^~'^du{u'^^^du) and the 0*-*''(m) have been given in the preceding subsection. 
We now make the following assumption, following the idea suggested in [28] : If we want 
to have scaling behaviour, we need that 

[ drt;(r) (4.30) 
Jo 
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at least for long times, which implies ^t,s ~ ~Ttx'^oit'^ — ■s'^)- Evaluating the expression 
flCTj) . we find 

se£" e,n=o 

^ r(2£ + i-2)r(2£+|-3 + f) 3 
r(2£ + f -2-n)r(2£+f -n-3 + f) ^ ' ^ 

Now we introduce u = \r\(t — s)~^^^ and .^t^^ ~ ~Tt\^o^^'^ ~~ ^'^^ ^^^^ expression. 
By changing the summation variables from n and i to k := 4i — 2n + s — 6 and i one 
can then read off the dynamical exponent, which is given by 

z = J. (4.32) 

The critical exponent z is therefore determined by the behaviour of the potential v{t). 
However, as we shall see later when treating the correlation function, only a value of 
F = \ leads to scaling behaviour of the correlation function. Nevertheless we keep F 
arbitrary for now and proceed with the autoresponse function R{t, s) = R{t, s; 0), which 
can be obtained from fl4.3ip by setting u = 0. All nonvanishing terms have to satisfy 
the condition 4i + s — 2n — 6 = 0, which excludes in particular odd values of s. Therefore 
only the contributions from the values s = 2 and s = 4 remain and the result for R{t, s) 
is 

Ro{t,s) = (t - s)-(^+i)/2 [d2g^'\t,s)+c,g^'\t,s] . (4.33) 
C2 and C4 are parameters and the expressions g^'^\t, s) and g^'^^t, s) are given by 

On- H--r^) 

g"\t, s) = (-4Ag.(( - s)-')hF, (I + J - |; |; (4-34) 



Note that we can write 
r2 



s''-'{^ ^ (4.35) 

t-s \ y-1 J 

with y = t/s. For future extension it is instructive to consider the different asymptotic 
behaviour implied by F, which we carry out in appendix C. Here we check only that our 
theory is in line with the exact results of the conserved spherical model as derived in 
section 2. If we choose C2 = ~7q^dTY\^ = Q-d = 0, we have seen in the last 

subsection that the expression Vr^F^'^\t, s, r) can be written as an integral (see fl4.27p ). 
On this integral representation, we apply formula (I4.28p . We also recall the fact that 
exp(^t^sV^)e^'^'' = e^^^-^t,!'^'^ _ It then follows from a straightforward computation that 
the result fl2.14p is reproduced correctly for kq = —Qa- 
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5. Response and correlation function in the noisy theory 

5.1. The response function 

We use the decomposition (14. ip and expand around the noise-free theory. Because of 
(14. lip , we find that the response function is equal to the noise- free response function, 
derived in the last section 

R{t,s;r) = Ro{t,s;r) (5.1) 

Therefore, we can take over the results already discussed in the last section. In 
particular, we can conclude that the form of the two-time response function in the 
critical spherical model with conserved order-parameter agrees with the prediction of 
local scale-invariance. 

5.2. The correlation function 

For the correlation function, the following terms remain (as usual r = x — y) 
C(t, s;r) = - ^j dudR {<p{t, x)<P{s, y)4>{u, R)V%^iu, i?))^ 

+ j^R {<Pit, ^)<Pi^, i/)0'(o, R))^ (5.2) 

where we have assumed uncorrelated initial conditions, that is 

{(f){R, 0)(j){R', 0)) = ao6{R - R!) (5.3) 

We denote the first term on the right-hand side of (15. 2p by Ci(t, s; r) and the second 
term by C2(t, s\r). Unfortunately, we do not have an expression for the three- and four- 
point functions. However, we can use Wick's theorem and the Bargman superselection 
rule, which leads to 

Ci(t, s- r) = / /^^ ('^('^' R))j^R{'t>{s, y)^{u, i?))^ 

^ du j dRV],(^4>{t,x)^{u,R))^(^^{s,y)^{u,R))^ (5.4) 



T 



16A 
and 



C2{t, s; r) = aoJ dR (^<p{t, a;)0(O, R)) J^<p{s, y^Q, i?))^ (5.5) 

Under the integrals, we always find two sorts of factors: The two-point function 
F^'^\t, s;r) = {(f)(t,x)(f){s,y)) was computed in appendix A and reads [r = x — y 
and u = \r\{t — s)^^^^) 

F^^\t, s; r) = {t- s)-"/' ^ Cs(j)^'\u) (5.6) 
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and the response function RQ(t,s,r) = {(f){t,x)(Vy)'^(p{s,y)) as computed in the 
previous section 

Ro{t,s;r) = (t - ^ c (5.7) 

When we introduce the expressions (15.61) and (15.71) with the appropriate arguments into 
(15 .4^ and (15. 5p we get the most general form of the correlation function fixed by LSI 
with z = 4. 



We now show, that this result is compatible with the exactly known result for 
C{t,s;r) of the conserved spherical model as derived in section 2. From the response 
function, we have already seen that we need to make the choice C2-d = c^-d = and 

C2 = — / j 3 %N y. Then we have the integral representation (14.271) for R(t,s;r) to 

"( 4 + 2 ~ A j 

which one can still apply the gauge transform (14.151) . For F^'^\t,s,r) one can find in 
the same way an integral representation, so that we have the following two expressions: 

dfc - 2-^ / ik ■ r 



Mt, s; r) ={t- sr^^^^y^ I ^ [ey-"^ exp 



(t-s)l/4 



Using these representations, we obtain 

Ci(t,s;r) = -— [\e{y - of^ + '^-^^Hl - of^ + '^-i'^ (5.10) 



2+22 I dk / ik ■ r \,j2\3~^ 



X (, + 1 - 2.)^.-- / (^ -P i - ^7773^ ) 



( ,4 t>' + sF-2{seY 

C2it,s;r)=aoit + sr-^''^^{ts)^'-^ ^ V (2^ """H " (^T^) 

f + s^ \ , , 

X exp + TT-r/tOT-- — ^ (5.11) 



16A "(t + s)i/2 

where we have used the fact that the integration over R gives a delta function. In order 
to compare with the exact results eqs. (12. 10112. lip , we see directly that we have the 
equalities Ai{t,s;r) = 6*2(^,5; r) and A2{t,s]r) = Ci(t, s;r) for x = |, A = — ^ and 
^ = I and r = \- Our symmetry-based approach has thus reproduced all the results 
from section 2, up to an identification of parameters. 
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6. Conclusion 

Current studies on non-equilibrium systems often concentrate on the scaling aspects, 
since then universality may be used to justify the analysis of rather artificial-looking 
models in order to obtain physical insight. In this work, we have been exploring the 
idea that dynamical scaling might be generalisable to a larger algebraic structure of 
local scale-transformations. This kind of prediction is most conveniently tested through 
the form of the two-time response functions and, since the form of the autoresponse 
function contains the dynamic exponent z only in combination with other exponents, 
the scaling form of the space-time response function will provide non-trivial information 
about how to construct local scale-transformation beyond the case z = 2 which by now 
is rather well understood. Since values of z quite distinct from two can be obtained 
for a conserved order-parameter, this motivated our choice to study the kinetics of such 
systems. The spherical model, quenched to T = Tc from a fully disordered initial state 
and the MuUins-Herring model with conserved noise are useful first tests, since their 
scaling behaviour is non-trivial, yet the models do remain analytically treatable. 

By looking for generalised local scale-transformation, it might appear at first sight 
that the non-invariance of a stochastic Langevin equation under such transformations 
would make such a programme futile. However, generalizing a technique from 
Schrodinger-invariance (which applies for example to phase-ordering kinetics with a 
non-conserved order-parameter ^8]), we have shown that also in the conserved case one 
can decompose the Langevin equation into a deterministic part which does admit larger 
symmetries and a noise part which breaks them. This decomposition is such that all 
interesting averages can be exactly reduced to averages calculable from the deterministic 
part of the theory. We have shown that this decomposition can be carried out in the 
critical spherical model and did confirm that the predictions obtained for the two-time 
response and correlation functions are fully compatible with the exactly known results 
in the spherical model. Together with the Mullins-Herring growth model, these are the 
first analytically solved examples which confirm LSI for a dynamical exponent z = 40 
In table 1 we collect the values of the exponents of the conserved spherical model and 
the conserved and the non-conserved Mullins-Herring models. For comparison we also 
list some of the values for the exponents of the bond-diluted 2D Ising-model, for which 
it was shown that the autoresponse function is in agreement with LSI [33]. 

That confirmation was possible because the deterministic part of the Langevin 
equation is still a linear equation. Numerical simulations in models where this is no 
longer so will inform us to what extent LSI with z ^ 1 can be confirmed in a more 
general context. At the same time, it will be necessary to derive a generalisation of the 
Bargman superselection rules in order to have a model-independent justification of the 
decomposition of the Langevin equation. Work along these lines is in progress. 

* Another example of a system with a dynamical exponent far from 2 which appears to be compatible 
with LSI is the bond-diluted 2D Ising-model quenched to T < Tc [33] . 
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model 


a 


h 




Ac 


Reference 


cons, spherical 


d/A-l/2 


d/A- 1/2 


d + 2 


d + 2 


|1] & this work 


cons. Mullins- Herring 


d/A - 1/2 


d/A - 1/2 


d + 2 


d + 2 


[1] & this work 


non-cons. Mullins-Herring 


d/A-l 


d/A- I 


d 


d 


liii 


diluted Ising 


0.24(2) 




1.32(4) 


1.280(20) 


[33] 



Table 1. Exponents a,b, Xu, of the conserved spherical model and the Mullins- 
Herring model at criticality and the bond-diluted 2D Ising model below the critical 
temperature. The latter model is defined by the hamiltonian Ti. — —J^iij) Jij'^i'^j 
with (Ti = ±1 and random variables Jij, which are equally distributed over the intervall 
[1 — e/2, 1 + e/2] with < e < 2. For the given values of the exponents of this model 
the temperature has the value T = 1.0 and the parameter e has the value e — 2.0, see 
|33j for details. For all models, the dynamical exponent is z = 4. 



Appendix A. The two-point function 



In this appendix, we compute the two-point function. We do the computation for d = 2 
for convenience, but the case for arbitrary dimension will be obvious. We write 



F^'\t,,t2,rt,rlr^„rl) = (^^(ti, r^, r^)02(t2, r^^ rf) 



(Al) 



where rf denotes the z-component of the coordinate vector and the y-component. 
Invariance under X_i, Y^l and Y^l gives us 



= F(t,r^r^), 
where we have defined: 



r^:=ri[-r2, := rf - r"^, t :-- 



Next we take the rotation generator 
and find 

{r'^drv -rJ'a,.)F(t,r^,r^) = 
From this we can conclude that 



tl-t2 



F = G{t, r) with r = 
Now we use invariance under Xq. It requires 



ri + r^y 



{tidt, + t^dt, + \ {rldr^^ + rfdrv + r^dr^^ + rld,{) + F = 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 
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with X = {xi + X2)/4. This can be transformed into 

tdt + ^rdr + x^G = (A8) 

Here we have used the facts 

= / o"^, o ^r, d^y = , ^ dr (A9) 



Before we proceed, we also note the fact that, when apphed to G{t,r), we have the 
identity 

+ dl = + d^y = dl + -dr = -dr{rdr) (AlO) 

'1 ^1 '2 ^2 f r 

This is just the operator in polar coordinates without angular part. For arbitrary 
dimension d we must use 

dl., +... + d% = -^^- drir^-'dr) (All) 

instead of (1A10I) . Turning to the operators Fi^''', we require 

4 

A = Ai = -A2, 7 = 7i = -72- (A12) 
Then invariance under Y^"'' eventually leads to 

ta,,.,) + 4Ar(->) (^^^dr{r^~^dr)^ - IQ^d,,.,, {^:^idr{r''-^ dr)^ j G = (A13) 

or, after multiplication with ^;p-y 

tdr + 4Ar ( -l-^dr{r^-^dr)^ - IQ^d, (^-^^^(^'^-i^^)^ j G = (A14) 

independently of the direction. Lastly, invariance under Xi yields after a straightforward 
but slightly lengthy calculation 



- 47ra, I G = (A15) 



Equations flASj) . flAlil) and flATSj) have to be solved. We multiply (jA8l) by -t and flATil) 



by — r/4 and add them to flAlSp . The resulting equation is satisfied provided that the 
condition 

xi = X2 (A16) 
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for the scaling dimensions holds. On the other hand (lASp is solved by 

G(t,r) =(5^,,^.,r(^i+^2)/V(M), with u = rt-^'^. (A17) 
From ( ]A14I) . the scaling function 0('u) satisfies the following equation 



This we rewrite using the shorthands a := 4A and h := — 167 

2 



5« (^^« (w'-'^n)) +au (^<9„ j + 6aJ = (A19) 



with 

-2 



u 

We note the following property 



:= ( [u^-'du) ) (A20) 



(^^« (^'"'^«)) w° = «(a + - 2)m"-' (A21) 

We make a similar ansatz as for the case d = 1 [11] 

00 

^c„m"+^ co^O (A22) 



00 

{u) 

n=0 



Introducing this into (lA19p . we obtain on the one hand, because of cq 7^ 



s G {0,2,4,2 -c/,4- 4 =: ^. (A23) 

On the other hand, we get a recursion relation for c„ 

_ a(n + s + c? - 6) + 6 ( \oa\ 

~ ~ (n + + s + - 2)(n + s - 2)(n + s + - 4)''""^ ^ ^ 

Third, we get an additional assumption on s, namely 

s ^ -2, -4, . . . , and s ^ -d, -d - 2, . . . (A25) 

If we did not have this, we would encounter an equation of the form c„ ■ + c„_4^ = 
with ^ 7^ 0, which would imply c„_4 = 0. However, this would be in contradiction to 
Co 7^ 0. The recursion flA24l) can be solved by standard methods (compare for instance 
[8]). One can show that for s G one may set ci = C2 = C3 = 00 and solve the recursion 
starting from Cq. We also set bi := c^e and obtain as the solution of the recursion: 



4; r(^-i + A) r(2£ + | + i)r(2£+f + f) 



(A26) 



jj This can in fact be deduced, except for some special cases in low dimensions. The latter are however 
covered by the different values of s 
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From this we can deduce the form of the scahng function (f){u) by equation (lA20p . The 
final result for the scaling function is therefore 

0(m) = ^c,0W(m) (A27) 

where Cg are free parameters and (p^^^u) are given by (substituting the original 
expressions for a and b) . 



oo 



0W(n) = 5^6f)M^^+^-4 (A28) 

with 



^ ^ ' r(2£ + f -i)r(2£+f + f-2) ^ ' 



This kind of function has been considered for instance in [311 132] and can in fact be 
rewritten as generalised hypergeometric functions with the result, up to normalisation 

, / Am^\ 2 ^ (d 731 d d Xu^ 

^"<") ~ (-1?) ■^'(4 + l-h'2n'4 + l^-Tr 

^'^"^^)-(-T^J ^^4'-V^-4'2-4'2^-l^ 

^ ^ ' V 16 y ^ V4 2 A' 2 4' 4' 2' 16 
Here the s-dependent prefactors have been dropped, as they can be absorbed into the 
parameters Cg. Note however that for certain values of s and 7/ A these factors might be 
zero so that the corresponding solution vanishes, see below for the free-field case. The 
coefficients as given in (IA291) will be used in the main text to determine the response 
function. Let us collect however some facts about the solutions (f)^^\u), as this might 
be relevant for models with non-conserving perturbations. 

• To be physically acceptable, the solutions have to be regular for m — > and van- 
ishing for u — >■ 00. The solutions belonging to s G {0,2,4} are indeed regular for 
M — > 0. However, the solution for s G {4 — d, 2 — rf} are not regular for d > A and 
have to be dropped for this case (for even d they where anyway excluded for d > A). 
For the same reason, the solution for s = 2 — d has to be abandoned for d > 2. 

Furthermore, we have bf^ = bf^^ and 6q°^ = 0, so that, up to a constant, the solution 
0'^^)(m) equals (f)^^\u). Therefore, we drop the solution belonging to s = 0. The 
general solution is then 

0(m) = ^c,0(^)(w) (A31) 
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with 

f {2,4} if c/>4 

£' ■= I {2,A,A-d} if 2<rf<4 (A32) 

[ {2,4,2 -c/,4- 4 if d<2 

For free-field theories, there is ^ = |. This makes the solution for s = 4 — d vanish 
(see flA29|) l 

The scaling function must vanish for u oo. Hence, not all of the coefficients Cs 
need to be independent, but they have to be arranged in such a way as to achieve 
this asymptotic behaviour. 



Appendix B. Calculation of an integral 

We outline briefly how to obtain the integral 



We choose spherical coordinates, that is: 
ki = kcos{9i) 
/c2 = sin(6'i) cos(^^2) 

; ; (B2) 

kd-i = ksm{6i) sin(6'2) ■ • • • ■ sin(6'd_2) cos(6'd_i) 
kd = ksm{9i) sin(6'2) ■ • • • ■ sin(6'd_2) sin(6'd_i) 

where k = \k\ and < < tt, . . . , < 6d-2 < ti", < 6d-i < 27i. Furthermore we rotate 
the system, so that we have k ■ r = kr cos{6i) with r = \r\. The Jacobian is given by 

• • • ' - k'-' sm'^-^e^) sm'-'{9,) • . . . ■ sin(^,_2) (B3) 



dik,e,,...,9d-i] 



If we change to the new coordinates in (IBll) . the angles 02, ■ ■ ■ , Od-i can be integrated out, 
yielding a prefactor A' , independent of r, a and (3, which we determine later. Therefore 

/•oo /"Tr 

j(-'/5)=^'/ dkk'^^^^^^e-'^^' dei(^''""''^''\m'''^{ei) (B4) 
Jo Jo 

The inner integral gives r(| — |)r(|)(ifcr/2)-'^~2/d_-^(iA;r), where Iv{z) is a 
modifled Bessel function pj], for which we use the series expansion Iv{z) = 
^^Q(n!r(z/ + n + l))~^{z/2y^'^'^. Then we perform the remaining integral over k. This 
leaves us with the expression 

^ n\T{l+n) \4 2 2 2) ^ ' 
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with a new factor A which is also independent of r, a and /3. This factor can now easily 
be determined by computing directly and then comparing with (IBSp . This yields 



A = 7i'^/^/2. Splitting the sum into even and odd terms and using for the Gamma 
functions the fact that r{2z) = 7r-2 2'^'-^r{z)r{z + i), we can rewrite this result in 
terms of hypergeometric functions and obtain 

Jk'* L 2 2'2'4'4 2'256ay 

8r(| + l + f) f \ fd ^ 3 d Id _r^ 

rfr (I + i + f) Vl6v^; ' ' U ^ 2' 2' 4 + 2' 4 + ^' 256a 
with the expression 



7r\ 2 r- _d_l_0 



27 ^ r(f + i)r(f) ^ ^ 

Appendix C. The response function in the scaling regime 

We reconsider equation fl4.33p and consider the different cases implied by F- 

• r = \ z = A . This case in fact corresponds to the conserved spherical model. In 
this case both g^'^\t, s) and g^^\t, s) depend only on the ratio t/s and the argument 
(I4.35I) of the hypergeometric functions tends to one for y ^ oo. Therefore, in the 
scaling regime we obtain: 

a = ^ - ^, and Xr = 2x + 2 (CI) 

• r < ^ z > 4: In this case, the arguments of the hypergeometric functions tend 
to zero in the scaling limit and we obtain 



Ro^,s) = s-^~^{y-l)-^^^ + ^^^^ 

if C2 7^ 0, this leads to 

i?o(t,s) = C2S-i-^(y-l)-f-5 (C3) 

and the same values (IClll for a and Xr as before. If C2 = 0, we have 

Ro{t, s) = d,s-^^+^-\y - l)-^-\y^ - 1) (C4) 

which means for the exponents: 

^ Xr X ^ 

«=2-r, -=2+l-r (C5) 

r > I < 4: If A < 0, the coefficients C2 and C4 have to be arranged in a way, 
that the exponentially growing parts of the hypergeometric functions cancel. More 
precisely, the condition 

Vil\ Til) 
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has to be satisfied. In this case we find 

R{t, s) = fos-f-^(^^-i)(i+t-i) [\^) (C7) 
where fo is a normahsation constant. This leads to the following critical exponents: 

. ^i-i + (2r-i)(i + ^-2) (C8) 

Ar/^ = i + ^ + (2/^ - 1) (l + i - (C9) 

If A > there is no condition on £2 and £4, as the exponential contribution decreases 
rapidly. In this case expressions flC7p and fICSp remain valid (only the constant fo 
may change) provided C2r(|) 7^ C4r(|). If the latter condition is not fulfilled, the 
result is 

(,r_i)x-(f-i) 



R(t, s) = fos-t-^(2r-i)(f-l) j v£ 1/ ] (CIO) 

(2/-I) 



and 



z 2 ' ' V4 A 

Ak/^ = ^ + ^ + (2r - 1) f 7 - ^ ) (C12) 



2 2' ' V4 A 

In all cases we find the relation 

Xji = z{a + 1) (C13) 



Appendix D. Conserved random walk 

Conceivably the most simple ageing system is the well-known Brownian particle, 
described by a Langevin equation only containing the noise and the external perturbing 
field [21] . Here we extend their idea to a conserved order-parameter and consider the 
following model, described by the Langevin equation 

dt<P{t,x) = -Vlh{t,x)+r]it,x) (Dl) 

where the rj is the same kind of noise (12.31) as in the main text and h is an external field. 
The solution of the model is promptly found and we obtain, where as usual r = x — y 

Rit,s;r)= -{Wl6ir))eit~s) (D2) 
C{t, s; r) = C(0, 0; r) - 2T(V^5(r))min(t, s) (D3) 

where the initial term C(0,0;r) can be dropped in the scaling limit. The limit 
fiuctuation-dissipation ratio is Xoo = 1/2. As in the non-conserved case [21], an 
equilibrium state is never reached. 
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The difference to the non-conserved case is that the space-dependence is described 
by a second derivative of a Delta-function instead of a simple delta-function. In order 
to better understand the meaning of this result, we return to a lattice discretization 
of the model. Then S{r) goes over to a Kronecker delta ^,,,0 and the second derivative 
becomes a discrete difference. We illustrate the result in one space dimension 



Then, for t > s, R{t,s;0) = 2 is a local maximum, while R{t,s;±l) = —1 would 
correspond to local minima. This is qualitatively similar to the form of R{t, s; r) seen 
in figure [2] for the conserved spherical model. 
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